Proton-Helium Spectral Anomaly as a Signature of Cosmic Ray Accelerator 
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The much-anticipated proof of cosmic ray (CR) acceleration in supernova remnants 
(SNR) must hinge on full consistency of acceleration theory with the observations; direct 
proof is impossible because of the orbit scrambling of CR particles. The recent ATIC, 
CREAM and PAMELA experiments indicated deviations between helium and proton CR 
spectra deemed inconsistent with the theory, since the latter does not differentiate between 
elements of ultrarelativistic rigidity. By considering an initial (injection-) phase of the diffu- 
sive shock acceleration (DSA), where elemental similarity does not apply, we demonstrate 
that the spectral difference is, in fact, a unique signature of the DSA. Collisionless plasma 
SNR shocks inject more He 2+ relative to protons when they are stronger and so produce 
harder helium spectra. The injection bias is due to Alfven waves driven by the more abun- 
dant protons, so the He 2+ ions are harder to trap by these waves because of the larger gyro- 
radii. By fitting the p/He ratio to the PAMELA data, we bolster the DSA-case for resolving 
the century-old mystery of CR origin. 



Cosmic rays (CR), discovered in 1912 [Q]], 
are subatomic charged particles with a power- 
law energy spectrum extended up to ~ 10 20 eV. 
At least to ~ 10 15 eV, they are commonly be- 
lieved to be accelerated by diffusive shock 
acceleration (DSA, or Fermi-I |2j-l6|]) mecha- 
nism, operating in supernova remnant (SNR) 
shocks (see lUS] for a review). Recent pre- 
cise measurements of proton and He 2+ spec- 
tra by PAMELA spacecraft [9J] indicate a small 
but significant difference between the two, 
confirmin g e arlier results of ATIC ifioll and 
CREAM lul E]. Since the DSA is electro- 
magnetic in nature and accelerates all ultrarela- 
tivistic species with equal rigidities alike, it was 
claimed inconsistent with this difference. 

Indeed, at the basic level the DSA mecha- 
nism predicts a power-law momentum distribu- 
tion oc p~i for the accelerated CR, where the 
index q depends on the shock Mach number 
q = 4/ (l — M~ 2 ). Therefore, q^ 4 — 4.1 seems 
to be rigorous for strong shocks (M ^> 1). At 
the same time, the subsequent escape from the 
Galaxy, partial escape of CR from the shock in 
the course of acceleration, and back-reaction of 
accelerated particles on the shock structure, in- 
troduce deviations of observed spectra from the 



above power law. Uncertainties in these correc- 
tions, not so much in the measurements, prevent 
validation of the DSA as the mechanism for the 
CR production in the Galaxy. 

Nevertheless, there is one fundamental prop- 
erty of this mechanism that can be tested inde- 
pendently from the above uncertainties. It is 
seen from the equations of particle motion in 
electric and magnetic fields E and B, written for 
the rigidity of CR nucleus 8% = pc/eZ, where p 
is the momentum and Z is the charge number: 
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Here &q — Am p c 2 /Ze, with A being the atomic 
number. These equations show that if protons 
and He 2+ ions enter the acceleration at 3% 3> 
Mq in a certain proportion N p /Nn e , this ratio 
is maintained in course of acceleration and the 
rigidity spectra are identical. Moreover, if the 
both species leave (escape) the accelerator and 
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propagate to the observer largely without col- 
lisions, they will maintain the same p/He ratio 
even if their individual spectra change consider- 
ably. 

The observations, however, were indicating 



for some time Jill, ll3 Ml6|l that the spectrum 
of He may be somewhat harder (by Aq = q p — 
#He ^ 0.1) than that of the protons over a wide 
range of rigidities Si ^> S%q. Recently, the 
PAMELA team ^ determined Aq = q p — qu e 
of Np/Nn e ratio as a function of rigidity with 
an unprecedented accuracy, Aq = 0. 101 ±0.001 
for Sg > 5 GV, where the finite Si effect fades 
out 113 311 . This finding challenges the DSA as 
a viable mechanism for galactic CR accelera- 
tion. The challenge is best seen from a remark- 
able similarity of the helium and proton spec- 
tra shown in an "enhanced" format, in which p- 
flux is multiplied by Si 2 - 8 and He-flux by Si 2J , 
Figffl 

While both spectra deviate from their power- 
laws, they do it synchronously (N p /Nn e is mea- 
sured with significantly higher precision than 
N p or Nu e , see below). First, let us focus on 
the following three common features of the He 
and proton spectra: (i) almost identical (three 
digits in the indices) convex shapes at 5 < M < 
230 — 240 GV with a likely roll-over towards 
the right end of this interval (ii) sharp dip at 
Si = 230 - 240 GV (iii) upturn with nearly the 
same slope at Si > 230 - 240 GV. 

These features are clues for possible accel- 
eration/propagation scenarios. In particular, the 
He and proton spectra cannot come from inde- 
pendent sources in their entireties. Otherwise, 
one is faced with the dip coincidence and the 
overall shape similarity. Neither can they come 
from a single shock, since the DSA and the sub- 
sequent propagation are inconsistent with the 
spectral variations shown in the features (i-iii). 
The remaining possibility seems to be that the 
low-energy part (M < 230 — 240 GV) originates 
from one source (SI) while the rest comes from 
the source(s) S2, including the invisible (un- 
der the SI) part with St < 230 - 240 GV. SI 
is likely to be a local source with soft spec- 
trum and a very low cut-off or a spectral break. 



The source(s) S2 generates a harder, feature- 
less spectrum that merges into (or comprises) 
the galactic background (see, however, Ill7ll for 
more scenarios). 

Despite considerable differences between 
the putative sources SI and S2, the p/He ra- 
tio is a remarkably featureless function of rigid- 
ity, oc Si \ in a wide rigidity range including 
the transition zone, at Si = 230 - 240 GV (see 
Ref-ISD and below). This points at a common 
(for S 1 and S2 and intrinsic to the DSA) mecha- 
nism that should account for the same 0. 1 dif- 
ference in independent sources. By virtue of 
eqs.£l]l2]), such a difference cannot arise in the 
region Si S?q. Therefore, it must originate at 
Si <C S$q, as we believe, in the following way. 

A small fraction of thermal upstream parti- 
cles, after crossing the shock may become sub- 
ject to the DSA (to be "injected") if they re- 
cross the shock in the upstream direction Ill8ll . 
Their amount depends on shock obliquity and 
Mach number (we will focus on quasiparallel 
shocks as more favorable for injection and fur- 
ther acceleration, but the results can be extended 
to the field inclinations w.r.t. the shock normal 
tf„ s ~30-40° HI ). 

In situ observations [20] of the Earth's bow- 
shock indicate that about 10~ 3 of incident pro- 
tons are injected. It is also known from such ob- 
servations that, on average, 1 .6 more He +2 ions 



than protons are injected [21]. This He/p injec- 
tion excess does not explain the PAMELA He/p 
excess unless it grows with the shock Mach 
number when the latter increases to the SNR- 
scales (M ~ 100). This is not known from in situ 
observations of shocks limited to Alfven Mach 
numbers Ma ~ Mjy 5. Therefore, we use the 



injection model II22I1 that predicts su ch g rowth. 
It is consistent with the observations 112 ill at low 
Mach numbers and with the recent simulations 
tl9|] in the important for the He hardening range 
of ~ 5 — 30. 

The mechanism of preferential He injection 
is based on the larger He gyroradius down- 
stream. Upon crossing the shock, both protons 
and He randomize their downstream frame ve- 
locity, which is ~ V s ( 1 — 1 jr) (where V s is the 
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shock velocity and r is its compression ratio) by 
interacting with magnetohydrodynamic waves, 
predominantly driven by the protons. We may 
consider the waves to be frozen into the flow 
since Ma = V s /Ca 3> 1, where Ca is the Alfven 
speed. As the proton gyroradius is a half of that 
of He 2+ (for the same velocity ~ V s ), the He- 
lium ions have better chances to return upstream 
since protons are retained by the downstream 
waves more efficiently. According to the model, 
the injection rates of both species decrease with 
Ma but the proton injection decreases faster. 

To quantify this effect, the model admits an 
initially unknown fraction of incident protons 
to return upstream where they drive a nearly 
monochromatic Alfven (magneto-sonic) wave. 
After being amplified by shock compression and 
convected further downstream, the wave traps 
most of the protons and regulates their return 
upstream. (He +2 ions are still regarded as test- 
particle minority). The monochromaticity of the 
wave upstream is justified by the narrowness of 
the escaping beam distribution compared to its 
bulk velocity upstream. The wave amplitude 
settles at a predictable level due to the obvious 
self -regulation of proton escape: if the escape is 
too strong, the wave grows to trap more protons. 

The mechanism is illustrated by Figfll where 
particle trajectories in the downstream wave are 
depicted in coordinates /i = Vji/V (cosine of the 
pitch angle w.r.t. the average magnetic field Bo) 
and a = k2Z + (j>, where k2 is the wave number 
downstream (related to that of the upstream by 
k2 ~ rk\), z is the coordinate (directed down- 
stream) parallel to Bo and shock normal, and 
is the gyrophase. Particles enter the downstream 
phase plane at its top when the shock sweeps 
in the negative a direction. Then they begin to 
move in the downstream wave along the lines of 
constant Hamiltonian 



the escape zone on the phase plane larger and 
more accessible to He 2+ ions than to protons. 
Note that in order to escape upstream, particles 
should cross the lines H = const which is en- 



abled by perturbations B4I1 

For Ma S> 1, the injection is suppressed ac- 
cording to rjp M A 1 In (Ma /A/*), where M* ~ 
10. The injection is more efficient at smaller 
Ma, but its dependence upon Ma is compli- 
cated [22]. A formal fit to the proton injection 
suppression factor gives r\ p 0.4 • M A p with 
o p ~ 0.6. The He 2+ injection is suppressed to 
lesser extent, yielding 7]He ~ 0.5 • M A aiie with 
OHe ~ 0.3. Both scalings are valid in the range 
5 ^ Ma i$ 100 so that the assumption about the 
test particle dynamics of He 2+ applies even for 
Ma ^ 100, where o p must be larger, according 
to the ri p oc M A l \n(M A /M*) asymptotic result. 
To accommodate this trend, we adjust o p within 
the range 0.6 - 0.9 and One within 0.15 - 0.3. 
Note, however, that very high Ma, where the in- 
dex o p grows, are linked with small SNR radii 
and their contribution is less important. 

From an SNR lifetime, we therefore select 
the Sedov-Taylor phase as the most important 
for the background CR production. The shock 
radius grows with time as R s ~ Csjt 2 ^, where 
C ST = (2.03£/p ) 1/5 , E is the SN energy and 
Po is the ambient density [23]. The shock speed 

is thus V s = (2/5)C 5 s ^R~ 3/2 . When the shock 
radius increases from R m i n to R mSLX , the follow- 
ing number of CRs (with momentum p) are de- 
posited in the shock interior 
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where B^ is the wave amplitude and v = 
kiV /(eZB/Ampc). For the same particle ve- 
locity V (which is an integral of motion), the 
parameter v for He, v He = 2v p , which makes 



where M is the current shock Mach number, 
M = V s /C s , a = p,He; C s is the speed of sound 
and the constant A is not important since we are 
interested only in the p/He ratio. The spectra 
can be represented as follows 



q(M) 
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Here 



-lnj 



is a reference (injection) rigidity, 



which can be arbitrarily fixed at M m ^ = 1 GV, 
since we are only concerned with the spectrum 
behavior at ffl > R inj , Re- 
introducing a new variable x = 
4\n(&/M m j), using the integration vari- 
able t — M~ 2 instead of M and substituting 
q = 4(l-M~ 2 ), r\ a oc M~ a ", for the p/He 
ratio we obtain 



N p /N Hs = C 



rb t o p /2 e -x/(l-t) dt 

j a 

f h t a He /2 e -x/(l-t) dt 
Ja 
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where the constant C is determined by the ratio 
of p/He concentrations. We also denoted a = 



M mlx ^ 1 and & 
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The result given by eq.© is shown in Figj3] 
along with the PAMELA p/He ratio. The agree- 
ment is very good besides the low rigidity range 
ffl < where it is not expected as the solar 
modulations, some further details of injection 
Il22ll . and possible but largely unknown propa- 
gation effects are not included in eq.©. There- 
fore, we make no attempts at fitting the & < 
~ ^inj range in Figj3l so the validity range 
of the fit, > R 2 , i.e., R > 2 - 3 GV, is 
clearly seen from the plot. The deviation from 
the highest rigidity point is likely to be due to 
large measurement errors and, in part, due to the 
breakdown of r\ a M~° a scalings. 
On representing Eq.© as 
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obtain for F 



4 In (R/R in j), we may 



Hp/No*** [\n{M/M^] ^ °^ 12 (7) 

The p/He ratio at ultrarelativistic rigidities, 
as opposed to the individual spectra, is not af- 
fected by the CR propagation, if collisions are 
negligible. Therefore, it should be examined 
for telltale signs intrinsic to the particle ac- 
celeration mechanism. The precise measure- 
ments of this ratio by the PAMELA [9], sug- 
gests reproducing their results theoretically with 
no free parameters. While we have obtained 
a convenient control parameter for this quan- 



tity, a(M A 



C7H e , from a collisionless 



shock model best suited to the PAMELA rigid- 
ity range, the model predictions need to be 
extended and improved systematically. Even 
though collisionless shocks is a difficult sub- 
ject of plasma physics, still not understood com- 
plet e ly [I24l -l26ll. we expect modern simulations 
II19l I27II to refine the proposed mechanism. 
This will extend the theory's fit to a broader 
range spectrum, currently being measured by 
the AMS-02, and help to determine whether or 
not galactic CRs are produced in SNRs. 

To conclude, there are alternative interpreta- 
tions of the He/p spectral hardening: (a) differ- 
ent SNR-type to contribute to the CR spectrum 
[0, S H9D, (b) variable He/p concentration in 
SNR environments 0, ED] and ic) CR spalla- 
tion 02I1 . They are reviewed in Ill7n . where it is 
pointed out that the overall data are best repro- 
duced if harder He spectra are directly released 
from accelerators. 
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where v = (7 a /2 + 1 and T denotes the gamma 
function. The last term in the braces, that cor- 
responds to the contribution from highest Mach 
numbers, may be neglected, as a <C 1 . For suf- 
ficiently large the p/He ratio behaves as the 
following power-law in In (M) 
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Figure 1: PAMELA fluxes of He and protons. He 
(solid line) and protons (dashed line) are multiplied 
by M 21 and M 2 /5, respectively (proton spectrum 
artificially reduced to emphasize its similarity with 
He spectrum). Circles represent PAMELA points 
adopted from Supporting Online Material for jjj. 
The sharp rise of the He beyond M ~ 800 GV is 
likely to be associated with growing errors, since 
it does not match with the ATIC-2 and CREAM 
JHEi] data at M > 10 3 GV. The "zig-zags" on each 
spectrum (also present at lower energies) are well 
within the error-bars (not shown here). 
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Fig.2 
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Figure 2: Protons (top panel) and He (bottom panel) 
in a monochromatic wave downstream. The trajec- 
tories (H = const, eq.0) are shown on the particle 
phase plane a, /J. for the same particle velocity V, 
yielding v = 1 .2 for protons and v = 2.4, for He. The 
wave amplitude B±/Bq = 4, which corresponds to 
the wave amplitude B± ~ Bq saturated upstream and 
compressed later by the shock. The vertical bars at 
the top of each panel schematically show the particle 
entrance from the shock surface when it moves to the 
left across the phase plane. 
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Figure 3: Main graph: PAMELA points (circles) 
and eq.© fit (line). The fit is obtained for a p = 
0.85, a He = 0.15, M min = 1.05, M max = 100 and 
the normalization constant C = 15.5 to match with 
the PAMELA data. This value is, however, con- 
sistent with the 0.1 He abundance. The twenty 
highest rigidity points are shown with the error-bars 
(stat.+syst.), where they seem to become significant 
and the rightmost point clearly deviates from the 
theoretical prediction (the data points adopted from 
the supporting online material of |@]). The proton 
and He spectral breaks (see also FigQ]), collocated 
(within uncertainties) at 230-240 GV, are shown with 
two vertical lines. At higher rigidities the data from 
ATIC-2 0] and CREAM H are shown in the in- 
set, however, as a function of energy per nucleon 
(both adopted from il2l1 ). 
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